Scaling Limit of Vicious Walkers, Gaussian Random Matrix Ensembles, 

and Dyson Brownian Motions 
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We study systems of interacting Brownian particles in one dimension constructed as the diffusion 
scaling limits of Fisher's vicious walk models. We define two types of nonintersecting Brownian 
motions, in which we impose no condition (resp. nonintersecting condition forever) in the future 
for the first-type (resp. second-type). It is shown that, when all particles start from the origin, 
their positions at time 1 in the first-type (resp. at time 1/2 in the second-type) process are identi- 
cally distributed with the eigenvalues of Gaussian orthogonal (resp. unitary) random matrices. The 
second-type process is described by the stochastic differential equations of the Dyson-type Brow- 
nian motions with repulsive two-body forces proportional to the inverse of distances. The present 
study demonstrates that the spatio-temporal coarse graining of random walk models with contact 
interactions can provide many-body systems with long-range interactions. 
FAGS numbers: 05.40.-a, 02.50.Ey, 05.50.+q 



The walks of independent random walkers, in which 
none of walkers have met others in a given time period, 
are called vicious walks. Since each random walk tends to 
a Brownian motion in the diffusion scaling limit, an inter- 
acting system of Brownian motions will be constructed 
as the scaling limit of vicious walkers Q, |^ . In an earlier 
paper, we explicitly performed the scaling limit of the 
vicious walkers in one dimension and derived the spatial- 
distribution function for the nonintersecting Brownian 
motions j^. 

The purpose of this Letter is to study the scaling limit 
of vicious walkers as a stochastic process, while the pre- 
vious paper reported its static properties at fixed times 
in order to clarify the relation with the timeless theory 
of Gaussian random matrices Q . We first claim that the 
limit process is in general time-inhomogeneous, i.e. the 
transition probability depends on the time interval T , in 
which the nonintersecting condition is imposed. Then 
we study in detail two types of processes in time t de- 
fined by setting T ^ t and T ^ oo, respectively. The 
former process is related with the Gaussian orthogonal 
ensemble (GOE) and the latter with the Gaussian uni- 
tary ensemble (GUE) in random matrix theory. This 
result demonstrates the fact that in the nonintersecting 
processes distributions at finite times depend on whether 
the nonintersecting condition will be also imposed in 
the future or not. We will derive the stochastic differ- 
ential equations for the latter type and show that the 
drift terms act as the repulsive two-body forces propor- 
tional to the inverse of distances between particles. In 
other words, the scaling limit of vicious walkers can real- 
ize the Dyson-type Brownian motion models (at the in- 
verse temperature f3 = 2) jsj. The Gaussian ensembles of 



random matrices can be regarded as the thermodynam- 
ical equilibrium of Coulomb gas system and that is the 
reason why Dyson introduced a one-dimensional model 
of interacting Brownian particles with (two-dimensional) 
Coulomb repulsive potentials. Here we should emphasize 
the fact that the vicious walkers, however, have only con- 
tact repulsive interactions to satisfy the nonintersecting 
condition. By taking the diffusion scaling limit, we can 
extract global effective interactions among walkers and 
they turn to be the long-ranged Coulomb-type interac- 
tions Q . Such emergence of long-range effects in macro- 
scopic scales from systems having only short-ranged mi- 
croscopic interactions is found only at critical points in 
thermodynamical equilibrium systems, but it is a typical 
phenomenon enjoyed by a variety of interacting particle 
systems in far from equilibrium. 

We also report the scaling limit of vicious walkers with 
a wall restriction. It will be shown that in this situation 
nonintersecting Brownian motions are discussed in rela- 
tion to the chiral (Laguerre-type) Gaussian ensembles of 
random matrices 0, [^i 

Let {^^^}fc>o, i G {lj2, •••,n}, be the n inde- 
pendent symmetric simple random walks on Z = 
{• • • , —2, —1, 0, 1, 2, • • •} started from n positions, 2si < 
2s2 < • • • < 2sn,Sj S Z. We fix the time interval as a 
positive even number m and impose the nonintersecting 
condition 



< < • • • < 1 < Vfc < m. (1) 



The subset of all possible random walks, which satisfy 
(|l|), is the vicious walks up to time m. For 2ei < 2e2 < 
••• < 2e„,ej G Z, let Vn{{Rl' }T=o; Rni = 2ej) be the 
probability that such vicious walks that the n walkers 



2 



arrive at the positions {2ej}^^i at time m are realized 
in 2"*" random walks. We then consider, with a large 
value L > 0, the rcscaled lattice Z/(-\/L/2), where the 
unit length is and let be the symmetric sim- 

ple random walk started from x on this rescaled lattice. 
The following was proved in the previous paper for 
given t > and < xi < X2 < ■ ■ ■ < XmUi < y2 < 
■■■ < Vn, liTaL^ooVniiRl'jkio'^RZ G , % + dy^]) = 
fn{t;{yj}\{xj})d'^y with 

fn{t; {%}|{x,}) = (2^i)~"/'^«to) (e^^/*, • • • , e^"/*) 
xe-E(-?+f?)/2* J] (e--^/* - e-./*), (2) 

l<j<fe<n 

where s\{zi, ■ ■ ■ , Zn) is the Schur function [b| and £,j{y) = 
2/n-j+i — (n — j) The function ^) shows the 

dependence on the initial positions of particles {xj} 
and the end-positions {yj} at time t of the probabil- 
ity that n independent one-dimensional Brownian par- 
ticles do not intersect up to time t. Its summa- 
tion over all end-positions is denoted by Nn{t]{xj}) — 
hi< - <v^ (^"'y fn{t\ {yj}\{xj}) and its asymptote for t 
oo was calculated as 



Afn{t;{x,}) ~ (2^)-"/2i-''"&„({x,})/c„ 



(3) 



with 77„ = n{n - l)/4, bn{{xj}) = s^(a.)(l, • • • , 1) = 
ni<j<fe<„(2;fc - Xj)/ik - j) and c„ 
ri!(23"/2]T'^^ r(j72 + where r(z) is the gamma 

function 

Since the vicious walkers are defined by imposing the 
nonintersecting condition (|l|) up to a given time m, the 
process depends on the choice of m. That is, the pro- 
cess is time-inhomogeneous. This feature is inherited in 
the process obtained in the diffusion scaling limit as ex- 
plained below. Let T > and consider the n noninter- 
secting Brownian motions in the time interval [0, T]. For 
< s < t < T, xi < • • • < a;„, j/i < • • • < ?/„, the transi- 
tion probability densities from the configuration {xj} at 
time s to {j/j} at t is given by 



]{s,{xjYt,{yj}) = 



fn{t^ s-{y,}\{x,})Nn{T -t-{y,}) 

Nn{T-S-{x,}) 



(4) 

since the numerator in RHS gives the nonintersecting 
probability for [0, T] specified with the configurations 
{xj} and {yj} at times s and t, respectively, and the de- 
nominator gives the probability only specified with {xj} 
at s, where we have used the Markov property of the 
process. The time-inhomogeneity is obvious, since RHS 
depends not only t — s but also T — s and T — t. 

In this Letter we study two special cases of T ; case 
(A) T = t and case (B) T ^ oo. We will show that there 
is an interesting correspondence between these choices of 
T and the Gaussian ensembles of random matrices Q . In 
order to see it we consider the limit Xj — > (1 < Vj' < n). 



Case (A) T ^ t. Since limt_o /„(i; = 
J^"^-!^ (5(a;j — yj) with Dirac's delta functions, 
Nn{0;{xj}) — 1 for any {xj}, and setting T — t 
makes (Q) depend only on t — s [Q. For i > and 



l<j<k<n 



X / d"j;s^(,)(l,...,l)e-E^?/2t (1 + 0(1^1)) 

'vi<-<Vn 

^n(n-l)/4 ^Xk/t _ pXj/t 



(27r)«/2c„ 



n 



l<j<k<n 



Then we have 
5^,(0, {0};i,{yj}) 



Cjit 



J 



{l + 0{\x\)). 



l<j<fe<ri 

< j/„ . If we set t — \ 



with C„ = n(n + l)/4 for j/i < 
and assume yi < ■ ■ ■ < ym then 

gi(0,{0};l,te}) = n! g^'''' {{y.}) , 

where the probability density of eigenval- 



ues of GOE random matrices |12 
Case (B) T ^ (X). Let 

Pn{s, {xj}; t, {yj}) = lim gl{s, {x-j}; t, {y^}). 

1 — >-oo 

By virtue of we can determine the explicit form for 
any initial configuration {xj} in this case as 

p„(0,{x,};i,{y,}) = (5) 

with the Vandermonde determinant hn{{xj}) = 
ni<j<fc<Ti(^fc ~ ^j)- particular, if we take the limit 
Xj ~^ (1 ^ Vj ^ n), we have 

l<j<k<n 

(6) 

with Cn = «V2 and ((27r)"/2 J]"^/ In this 

case we set t — 1/2 and assume yi < • • • < y„. Then 
p„(0,{0};l/2,{y,}) = n!.g^rUE(|y^|)^ 

where On^^ {{yi}) is the probability density of eigenval- 
ues of GUE random matrices. In this case the noninter- 
secting condition will be imposed forever (T ~* oo), while 
in the case (A) there will be no condition in the future. 
The distributions of particles at finite times depend on 
the condition in the future. 

By generalizing the calculation, which we did in the 
case (A), for arbitrary T and comparing the result with 
(^, we have 

g^(0,{0};t,fa}) 
p„(0,{0};t,{y,}) 
= c r-"("-i)/4^"("-i)/2 -^niT - t; {yj}) 



Ui<j<k<niyk-yj) 



3 



with c„ = c„/< = W22-"n;=ir(j + i)/r(j/2 + 1). 

This equahty can be regarded as the multi-variable gen- 
eralization of Imhof's relation between the probability 
distributions of Brownian meander and the Bessel pro- 
cess Jl3[ . 

In the case (B) we have obtained the explicit expres- 
sion (|^) for arbitrary {xj} and we can derive a system 
of stochastic differential equations for the process. Using 
it we will explain why we have the GUE distribution at 
time t = 1/2. Let 



1 



for A: = 1,2, 



• n. 



It is easy to verify that 

'^''{{xj}) = dklogh„{{xj}), 

71 

Y.[dka'{{x,}) + {a'^{{x,})y 



k=l 



0, 



(7) 

(8) 



where dk = d/dxk- Using these equalities, we will prove 
shortly that pn{0, {xj};t, {yj}) solves the diffusion equa- 
tion 



1 " 
'(i; {xj}) = {xj}) + a''{{xj})dkuit; {xj}), 



d_ 

dt 

(9) 

where A = X]fc=i^fc- This implies that the process de- 
fined in the case (B) is the system of n particles with 
positions X2{t), • • • , Xn{t) at time t on the real axis, 

whose time evolution is governed by the stochastic dif- 
ferential equations 

dxk{t)=a^{{xj{t)})dt + dBk{t), k = l,2,---,n, (10) 

where {Bk{t)}^^i are the independent standard Brow- 
nian motions; Bj{0) = 0,(Bj{t)) = and {(Bjit) - 
Bj{s)){Bk{t) — Bk{s))) = \t — s\Sjk- Because of the 
scaling property of Brownian motion, ^/KBj{t) is equal 
to Bj{Kt) in distribution. Then, if we set t — 2t' 
and write Xk{t) — Xk{t'), ( |l0|) is equivalent with 
the a = 0,(3 = 2 case of the equations dxk{t') = 
-l3dkW"'[{xj{t')})dt' + y/2dBkit'),k = 1,2, with 
W^{{x,}) = aE;=i^|/2 - Ei<,<A<„log(^fc - 2:,). 
When a = 1, they are known as the stochastic differ- 
ential equations for the Dyson Brownian motions at the 
inverse temperature (3 — 1/kT and the stationary distri- 
bution oc exp(-/3Tyi({xj}) |, |. If a = 0, the factor 



exp(— a;|/2) is replaced by exp(— ^ i|/4t') for finite 
t' and thus when t' = 1/(2/3) we may have the Gaussian 
distribution with (3. Setting (3 = 2 gives t = 2t' = 1/2. 
It should be noted that the system of diffusion equations 
describing the Dyson Brownian motions with (3 — 2 can 
be mapped to the free fermion model p^ . 

Now we prove that (^) satisfies (^). First we remark 
that 01 



fn(t;{yj}\{xj}) 



det 

l<j,fe<n 



(2TTt) 



-\I2 -(Xk-Vjf I2t 



That is, /„ is a finite summation of the products of Gaus- 
sian kernels and thus it satisfies the diffusion equation 
Therefore 



d 



g^p4^;{y.}|{a:,}) = 2/,„({,^.}) 



—^Mt; {yj}\{xj})- 



Then we can find that, if {a'^{{xj})} satisfy the equations 

n 

fc=l 



:dkfn{t;{yj}\{xj}) 



hni{Xj}) 

n 

J2 dk{l/hr.{{x,}))dkfn{t; {yM^A) (11) 



and 



n -. n 

Y^''({^j})Myh,,{{x,})) = --Ydl{l/hn{{x,})), 



k=l 



k=l 



(12) 

the proof will be completed. It is easy to see that (O) 
is satisfied if (0) holds for all k. Moreover, using 
we can reduce (|12[) to (^. Then the proof is completed. 
This proof shows that the origin of the long-ranged in- 
teractions (|^) is the normalization factor l/hn{{xj}) in 

(!)• 

In order to derive the chiral versions of the Gaussian 
ensembles of random matrices @, § and of the Dyson 
Brownian motions, next we consider the vicious walker 
problem with a wall restriction ||l6| . We impose the con- 
dition 



Rl' >0 1 < Vfc < 



(13) 



in addition to (|l]). That is, there assumed to be a wall 
at the origin and all walkers can walk only in the region 
a; > 0. We write the probability to realize the vicious 
walks which satisfy these conditions (|l]), (^3|) and {Rm — 
2ej}"^i as V+{{Rl'}'^^g;R'rA = 2ej). By the Lindstrom- 
Gessel-Viennot theorem and the reflection principle of 
random walks [061 , we have 



K^+(Knr=o;^^=2e,) = 2 
X det 



TO 

Sk - e. 



l<j,fc<n \ \ "2" 



Y+Sk + {ej + 1) 



Following the same calculation as was done in we 
will obtain, for < > 0, < xi < • • • < Xn, < yi < 
■■■ <yn, limi^oo V+i{Rl'ma;Rl] e [y,,yj + dy,]) = 
/,t(t;{%}|{3^j}K2/, with 



3 = i 



n {(' 



.Vk/t 



1)}, 



l<j<k<n 



4 



where sp\{zi,- ■ ■ , Zn) is the character of the irreducible 
representation specified by a partition A = (Ai, • • • , A„) 
of the symplectic Lie algebra (see, for example, Lec- 
tures 6 and 24 in @) given by det(zj''=+"^''+i - 

z-(^^-+"-"+'Vdet(z;-'^-+i - z7^""'=+')), and we set 

- i^ti^), ■ ■ ■ with = x„-,+i - (n - 



j + 1). It is known that spa(1, • • • , 1) = l\i<j<k<n 
^j)/(™fc~2Z?i) HjLi ^j/'^i with^j = Xj+n-j + l,,,oj 
and 



m, = 



n 



,2|27 



\<j<k<n 



n 



.+7«(n-i) -Q r(i + j7)r(a + 70--i)) 



r(i + 7) 



for Rea > ((17.6.6) on page 354 in g). Then 
we can obtain the asymptote of A/jf (i; {xj}) — 

/o<yi<...<y„ ^"y/n (<; {Vi)\{^i)) for i » 1 aS 

AA+(i; {X,}) = i-"(2»+i)/2(i + 0(1/,)) M-lll 

/n 



l<j</c<n 



2«(n+2)/2 ^^/^ j\ / j_l 



with ?7, 



+ _ 



72, where 6+({x,}) = Jli 



(4 



Two types of diffusion processes can be defined as in 
the cases without a wall. In the case (A) we will con- 
clude that at time t = 1 the positions of the nonintersect- 
ing Brownian particles all started from the origin with a 
wall restriction is identically distributed with the positive 
eigenvalues of random matrices in the chiral GOE. In the 
case (B) the transition probability density is given as 



p+(0,{x,};t,{yj}) 



h+i{yj}) 
hti{x,}) 



/n+(^;{%}|{2;,}) 



<j<k< 



We 



with h+{{x,}) = 111 

can prove that p^{0, {xj};t, {yj}) satisfies the diffusion 
equation (^ if a''{{xj}) is replaced by a'^+({xj}) = 
jyU^j^k ^^k/ixl - x]) + l/xk. When {x,} = {0}, the 
distribution at time i = 1/2 may be equal to the distri- 
bution of the positive eigenvalues of random matrices in 
the chiral GUE. 

In summary, we have studied the diffusion processes 
derived as the scaling limits of vicious walkers in one 
dimension with and without a wall. Two types of in- 
teracting Brownian particles were defined depending on 
the situations in the future and interesting correspon- 
dence between these processes and the Gaussian ensem- 
bles of random matrices was discussed. The systems of 



stochastic differential equations for the second-type pro- 
cesses were determined, which can be identified with the 
Dyson Brownian motion models at /3 = 2 with appro- 
priate change of time-scales. Further study of general 
time-inhomogeneous processes determined by (^) will be 
an interesting future problem. 
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